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Abstract 

We report on the study of extremely high energy photons, pions and 
etas, produced in intermediate energy heavy-ion collisions. Possibility 
of imaging the final-state phase space in these collisions by the Bose- 
Einstein correlations for photons is critically examined. 
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1 Introduction 



Measurement in heavy-ion (HI) collisions of hard photons of energies several 
times higher than the beam energy per nucleon, has stimulated large interest 
in the development of an appropriate theoretical framework for the descrip- 
tion of rare energy fluctuations which could account for the measured proper- 
ties of particles. In this context, the success of the stochastic multiplicative 
processes in describing the subthreshold particle yields came as a surprise. 
Nucleon in the transient state of the HI collision is described by the stochas- 
tic equation : dp/dt = F + G(p)X , where A is a random Gaussian process 
(< A >= 0, < A 2 >= a 2 ) accounting for fluctuating nuclear environment, and 
G{p) = 9o + 9iP + ' " ' is a smooth function of the particle momentum. In the 
leading order, the time evolution of the relative momentum of two particles is 
: d(5p)/dt ~ a(5p)\ and, hence, the relative kinetic energy evolves as : 

-^Enn ~ 2aE NN \ . (1) 

The relative energy of the nucleon-nucleon pair in Eq. (jl|) is a stochastic, 
short-correlation quantity driven by the multiplicative noise [l], 0] . Taking e(t = 
0) = 5(Enn — Ecm/A) as the initial condition, one may derive the limiting 
probability distribution of the relative energy field : 

^-(-te)!) , 

where N is the normalization factor, yfs is the width of the Gaussian variable 
and /i = \og(EcM /2A) + s/2 . The production cross section of a particle V 
in the NN collisions is then : 

* (a ° = r*deP{e)aP N {e-E$) , (3) 

th 

(x) 

where E\ h is the energy threshold in the NN collision and E tot is the total 
energy available for the production of a particle in the HI collision. Assuming 
that the time-duration of instabilities is approximately independent of bom- 
barding energy in a broad energy interval and that the parameters of random, 
short-correlated medium depend weakly both on energy and on the kind of 
produced particle, it was shown that all data on total and differential cross 
sections of produced pions at Ecu < lOOMeV/A is remarkably well fitted by 
a distribution (H) with 2 parameters : iV and s in the whole region Jl|. Also 
the spectrum of deeply subthreshold 77's at E C m/A = 180 MeV is reproduced. 



The above stochastic multiplicative process is an attempt to characterize 
the global features of relative energy fluctuations in the off-equilibrium media 
without entering into details of elementary processes involved in the produc- 
tion of particle x. Alternatively, one can start from the cascade picture of 
particle-particle collisions, including detailed informations about all elemen- 
tary processes in the medium and, in this way, one can try to reconstruct the 
particle production scenario in HI reactions. Below, we shall use for this pur- 
pose the Dubna Cascade Model (DCM) which has been recently extended 
by including various radiative capture processes such as : np — > d^, ir, r\ ■ ■ •) , 
as well as the off-shell evolution of nucleons and A'sjH . 



2 Subthreshold pion dynamics as a source of hard pho- 
tons, pions and etas 

We shall omit discussion of ingredients concerning elementary processes con- 
sidered in DCM. Interesting reader can find further informations in Refs. [[|. || . 
In the calculations of the photon spectrum, we consider the pn bremsstrahlung 
process in the one-boson approximation and we omit the pp bremsstrahlung. 
We include also the radiative capture process : pn — > d'-f , which can be impor- 
tant in the description of high energy part of the photon spectrum. This pro- 
cess requires the existence of bound state of final deuteron, i.e., the deuteron 
momentum must be above the Mott momentum if the influence of the sur- 
rounding nucleons is taken into account^ || . The process pn —>■ d'-f is most 
important at E/A ~ 150MeV[|], |S| . The decay processes : 7r° — ► 77 and 
A — > N'-f , which could give rise to hard photons, have been included as well. 

It turns out that the most important source of hard photons is the process : 
ttN — * iV7 ||. Pions in the DCM are produced either directly : NN NNir , 
or in two steps involving intermediate A or off-shell nucleon|| : NN — > AN 
or NN -> iV (o//) iV and then A — > Ntt , AiV — > NNir or N (of ^N -> A^iVvr . 
Formation of A is calculated so that the effective mass of the irN system follows 
the A mass distribution, its width being dependent on the pion momentum in 
the 7rN system. The yield of the primordial pions is subsequently modified by 
the absorption and rescattering processes in the nuclear medium[||]. The same 
two-step mechanism as well as the process ttN^ ^^ — > Nn^n- ■ •) have been 
considered as a source of hard pions and deeply subthreshold etas. 

Fig. ^] shows the comparison of the experimental data for the 7r + production 



in p + Ar reaction [jnj in the broad range of energies with the results of DCM 
calculations. The agreement between theory and experiment becomes even 
more impressive if one notices that the experimental data contains a systematic 



error of about 40 % [HTO 



Fig. [T] compares calculated and measured photon spectrum (full points) 
in the reaction 181 Ta + 197 Au at 39MeV/A || The dominant contributions 
: pn — > pwy , pn — » dj and Nrr — > Nj , are shown for the comparison. 
The contribution from the decay 7r° — > 77 is insignificant at this bombarding 
energy. The agreement between DCM calculations of the photon spectrum 
and the experimental data is excellent up to the highest photon energies. 

Finally, Fig. [1] presents measured |TT[| (full triangles) and calculated (full 
points) transverse mass spectrum of pions for the reaction Ar+Ca at 180MeV/A. 
One may notice the deficit of hard pions (m t > 0.4GeV ) in DCM calcula- 
tions. The measured spectrum decreases slower than exponentially, similarly 
as in the log-normal behaviour. The open points in Fig. [I] show the contribu- 
tion from processes involving off-shell nucleons. The crosses give the calculated 
^-production cross section. The estimated meson ratio a^/a^o is < 3 • 10~ 6 in 
the data as compared to 3 • 1CT 7 in the DCM[|TT| . 



Imaging the final-state phase space in HI collisions by 
HBT correlations 



Here we want to examine the status of HBT measurements^] as a tool to 
study the spatio-temporal extension of the emitting source at intermediate 
energy HI collisions. The two-particle correlation function is defined as : 

C ex (r) \ = A(pi,p 2 ) = I Jexp(iqx)S(x,K)d 4 x\ 2 

2 {PUP2) ~ Pi(pi)PM [JSixup^x,] ■ [JS(x 2 ,p 2 )d*x 2 ] ' 1 1 

where S is the source term (the Wigner transform of the single particle density 
matrix), q = p 2 —p\ and K = (1/2)(E 1 + E 2 , p x + P2). Cfj* is positive definite 
even though S(x,p) may take both positive and negative values. Eq. ([|) 
constitutes a formal basis for incorporating the Bose-Einstein symmetrization 
in the cascade codes of HI reactions. However, in the actual applications, one 
replaces this exact result bypX 



C app (p p ) - I / exp(iqx 1 )S(x 1 ,p 1 )d 4 x 1 J exp(iqx 2 )S(x 2 ,p 2 )d 4 x 2 | 
2 ' [J S(x 1 ,p 1 )d 4 x 1 ] ■ [J S(x 2 ,p 2 )d 4 x 2 ] 

which corresponds to weighting the pair at space-time points (xi,pi) , (x 2 , p 2 ) with 
l + cos[(pi— p 2 ){xi— x 2 )] , and is equivalent to the assumption of smoothness of 
the source function. This assumption may lead to certain pathological results 
(see, e.g., the discussion in Ref. WW) ■ 



The advantage of photons over other probes is that they escape from the 
interaction region unperturbed. Let us take the photon source in the form : 

S(x, K) = Nexp{-E/T) exp(-r 2 /i? 2 ) exp(-t 2 /r 2 ) , (6) 

which approximates well the source function extracted from the DCM evolution 
[f|] . The correlation function corresponding to (§) is : 



(7) 



where z = Q/T , Q 2 = q 2 — g 2 and 1^2(2) is the McDonald function. The 
polarization factor JT^| : A = (1/4) (1 + cos 2 Oi 2 ) , where 6i 2 is the angle 
between photon momenta, has been taken equal 1/2. Using (|5|) one finds : 

C? P (Q, go) = 1 + A exp ( - ^) exp ( - |(r 2 + i? 2 )) . (8) 

Note the absence of the pre-exponential function, which is a measure of im- 
portance of the off-shell effects and disappears if : Q — > or T — > oo. 

Both exact (|7|) and approximate (H) correlation functions, are plotted in 
Fig. |H for two source sizes : i?o = 4fm and 8fm, and for the 'apparent temper- 
ature' T = 15MeV, which is a typical value in HI reactions at E/A < lOOMeV. 

for i?o = 4fm shows a strong increase due to the fact that a typical photon 
wavelength in this low-T regime is comparable to the source size. For both val- 
ues of Rq the approximate correlation function provides a poor approximation 
of the exact one, i.e., the cascade calculations using the smoothness approx- 
imation (j5|) are not reliable. These calculations could be improved using the 
smearing prescription (T^] where the 'classical' source is convoluted with a min- 



imum uncertainty Gaussian wave-packet : g(x,p) ~ exp[— x 2 /(2cr 2 ) +2p 2 a 2 )} . 
However, in this case results both for Ci and for the energy spectra depend 
on the value of the smearing parameter o 2 . 

In the experimental analysis, the source size could be extracted more reli- 
ably using Eq. (0) in the fitting procedure. One should however notice that 
for small source sizes, the rise of the correlation function can be easily con- 
fused with a Gaussian peak at Q ~ m n from the decay : 7r° — > 77 . In this 
unfavorable case, the experimental extraction of Rq may be very difficult. 



4 Conclusions 



We have shown that cascade dynamics with the schematic account of mean- 
field dynamics, is successful in reproducing spectra of extremely high energy 



photons and hard pions produced in the deeply subthreshold regime of HI 
reactions. The significant deviations between the DCM calculations and the 
data begins to be visible only at extremely large m t 's. More data in this 
domain, which may provide a stringent test of the time-evolution of energy 
fluctuations in the off-equilibrium nuclear medium, are badly needed. 

In principle, the same cascade scenario could be tested also in HBT mea- 
surements. However, the rapid shape variations of the 77- correlation function 
with the source radius and the importance of off-shellness (Q 7^ 0) in the in- 
termediate energy HI reactions may pose serious problems in extracting the 
quantitative information. HBT measurements using pions in this domain are 
also not free from these defects but the corrections in this case are somewhat 
easier to introduce both in the calculations and in the experimental analysis^] . 
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Fig. 1. Experimental energy and angle integrated yield of tt + from p + Ar reactions 
is compared with the results of DCM (stars). 
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Fig. 2. Measured photon spectrum in the reaction 181 Ta + 197 Au at 39 MeV/A 
is compared with the DCM calculations (the solid line). The calculated spectrum 
is decomposed into fractions corresponding to the following elementary processes : 
pn — > pnj, Nir — > TVy, — > tZ-y. 
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Fig. 3. Experimental pion transverse mass distribution in the reaction Ar + Ca at 
180 MeV/A is compared with those obtained in the DCM. For details see the de- 
scription in the text. 



S(E) ~ exp(-E/T),T=15MeV 
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Fig. 4. Four different projections of the correlation function calculated analytically 
for the source distribution (||). The solid lines show results obtained exactly from 
the definition (|j) for the two source sizes : Rq = 4 and 8fm, whereas the dashed and 
dotted curves have been obtained in the smoothness approximation (||) for Rq = 4 
and 8fm respectively. 



